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The mean field effect in quark-gluon plasma on J/il) production in heavy ion collisions is pertur- 
batively calculated in a transport approach. While the global nuclear modification factor Raa is 
not sensitive to the mean field, the reduced threshold for J/tp regeneration leads to a significant 
enhancement of Raa at low transverse momentum at RHIC and LHC energies. 
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J/ip is a tightly bound state of heavy quarks c and c, its dissociation temperature Td in hot medium is much higher 
than the critical temperature Tc of the deconfinement phase transition. Therefore, the measured J/ips in nuclear 
collisions at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) carry the information of 
the hot medium created in the early stage of the collisions and has long been considered as a signature 1] of the 
new state of matter, the so-called quark-gluon plasma (QGP) 2]. 

The hot nuclear matter effect on J/tp production includes elastic and inelastic interactions. The latter, for 
instance the gluon dissociation process g + J/ip c + c and its inverse process, leads to a J/ip suppression |3l4ll| 
and regeneration [l3 - [l7| and is widely considered in the models for J/ip production in heavy ion collisions. The 
former, which changes the particle properties like mass, is weak for a heavy quark system at low temperature 
and therefore neglected in almost all the models. However, from the calculations with quantum chromodynamics 
(QCD) sum rules [11] and QCD second-order Stark effect there is a sudden change in the mass of J/tp around 
Tc- For instance, at temperature T/Tc = 1.05 the mass shift 5mj/^ — mj/^{T) — mj/^{Q) can reach —(100 — 200) 
MeV, which is already comparable with the mass change for light hadrons and should have remarkable consequence 
in Jf-ip production in heavy ion collisions at RHIC and LHC where the fireball temperature is believed to be much 
higher than Tc- 

Let's first qualitatively estimate the mean field effect on the J/tp distribution. The decreased mass reduces the 
threshold value for the J/tp production in QGP, and should result in an enhancement for the J/ tp yield. Secondly, 
the attractive force between the inhomogeneous medium and J/tp, 

F(x, p) = -^Vm,/, = _!!^^vT (1) 



with J/tp energy Ej/^ — y "^j/^ + P^j pulls J/tp to the center of the fireball and kicks J/tp to a lower momentum 

region. This will lead to an enhancement at low momentum and a reduction at high momentum. In comparison 
with the reduced threshold by the temperature dependent mass mj/^(T), the mean field force induced by the 
inhomogeneous fireball T(x) is a second order correction. For a homogeneous fireball, the threshold is reduced, 
but the force disappears. 

Since charmonia are so heavy and difficult to be thermalized in hot medium, we use a detailed transport 
approach to describe their distribution functions /^(p,x, rjb) in the phase space at fixed impact parameter 
b of a nucleus-nucleus collision, 

^ + ^ • V,/* + F • \7ph = -a^h + (2) 

where V^; and Vp indicate three dimensional derivatives in coordinate and momentum spaces. Considering the 
fact that in proton- proton collisions the contribution from the decay of the excited charmonium states to the J / tp 
yield is about 40% [2l|, we must take all the transport equations for ^ = J/tp,tp' and Xc in the study of J/tp. 
The collision terms on the right hand side are from the inelastic interaction with the medium. The lose term 
a [13, [131 arising from the gluon dissociation process controls the suppression, and the gain term /3 [13] related to 
a by detailed balance governs the regeneration. The cross section for the gluon dissociation axs,{0) in vacuum has 
been calculated through the method of operator production expansion 22, 23] , and its value at finite temperature 
can be obtained from its classical relation to the charmonium size r<f, (t*(T) = (r,T,(T))/(r,T,(0))g^(0), where 
the averaged r-square can be derived from the Schrodinger equation with lattice simulated potential j25l, [Hj . 
Different from the previous study, we consider here not only the inelastic processes, but also a mean field term 
F • Vp/*. The local temperature r(x, r) appeared in the elastic and inelastic terms are obtained from the 
hydrodynamic equations [10| for the space-time evolution of the QGP. The charm quark mass rric is determined by 
the charmonium mass m^j and its binding e nerg y [3], rric = {mqi — e*) /2. The cold nuclear matter effects on 
J/tp production, such as nuclear absorption '27], Cronin effect [H, [2^ and shadowing effect can be reflected 
in the initial charmonium distributions /if (p, x, To]b) with tq being the formation time of the hot medium. To 
simplify the numerical calculations, we neglect in the following the nuclear absorption and shadowing effect at 
RHIC and LHC energy and consider only the Cronin effect with a Gaussian smearing treatment [l3, [lH . 
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The transport equation with the mean field force can not be solved analytically. Supposing that the very small 
elastic cross section for charmonia at hadron level is still true at parton level, the mean field term can be considered 
as a perturbation, and the transport equation can be solved perturbatively. Fully neglecting the mass shift, the 
zcroth-order transport equation can be strictly solved with the solution [20| 



/4°Hp,x,r|b) = /*(p,xo,ro|b)e 



- dTiQs,(p,Xi,Ti|b) 



+ 



/ dri/3*(p,xi,Ti|b)e •'^i 

J Tn 



(3) 



with the coordinate shift x„ = x — p/£'^(t — t„) which reflects the leakage effect [l|, |32] during the time period 
T — r„ . The first term on the right hand side is the contribution from the initial production, it suffers from the 
gluon dissociation during the whole evolution of QGP. The second term arises from the regeneration and suffers 
also the suppression from the production time ti to the end of QGP. 

Taking into account the fact that the dissociation temperatures for the excited states are around Tc 2J|, we 
do not consider their mass shifts in QGP. Therefore, for the transport equations for ijj' and Xc, their masses 
are temperature independent and the mean field forces disappear, the zeroth-order solutions become the full 
distributions, = and /^^ = f^J. 

We now consider the mean field effect on the J/ip distribution fj/^. We extract the mass shift from the 
calculation using QCD second-order Stark effect jl9| . 



5mji^{T) = mji^,{T)~mji^{Q) 



' 18 e 



(4) 
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with constants a = 0.271 fm and e = 640 MeV, the temperature dependent coupling constant a^, and the functions 
Mq = e — 3P and M2 = e + P determined by the energy density e and pressure P of the medium and extracted 
from lattice QCD simulations js^]. Leaving out first the mean field force induced by the inhomogeneous property 
of the fireball but keeping the temperature dependence of the mass, the solution of the transport equation has the 
same form with the zeroth-order distribution, the only difference is the replacement of toj/^(0) — >■ mj/^{T), 



fW ^ AO) I 



(5) 



If the effect of the mean field force is not very strong, we may infer that it is not such a bad approximation to 



replace /j/^ in the mean field term by the known /j^^, the approximate transport equation 

dfj/^ p 



^xfj/^ — —oij/^fj/^ + Pj/^ — F • Vp/j^^ 



dT 



E 



is then similar to the equation for /i/L but with a replacement for the regeneration 



(1) 



(6) 



(7) 



This means that the mean field force can be considered as an effective regeneration, which is not from the 
coalescence of heavy quarks but due to the J/ip mass shift. The approximate solution so obtained is known 
as the second-order J/'ip distribution. 



/(5;(p,x,r|b)=/W(p,x,r|b 



dri 



F(xi).Vp/,W (p,xi,ri|b) 



— J (ir2Q(p,X2 ,r2 |b) 



(8) 



(2) 

Substituting the obtained fj/^ into the mean field term, and solving the transport equation again, we can 

derive the third-order distribution function f j^[p- With the similar way, the distribution to the n-th order can be 
generally expressed as a series of the mean field force, 



n-l 



/r(p,x,r|b) = /4^)(p,x,T|b)- ^(-ir 



dro - ■ ■ 



dT„ 



nF(x,).vJ /4')(p,x„„T„,|b) 



— r dr' a{p,x' ,t' \h) 



(9) 



Note again that the mass shift affects the J/ip production in two aspects. The reduced threshold, which is 
considered already in the first-order distribution /j/l^, makes the production more easily, and the attractive force 
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between J/ip and the matter, which is included only in the higher-order distributions /^"^ 
a momentum shift and leads to a low pt enhancement and a corresponding high pt suppression. 

By integrating the charmonium distribution function /.(.(p, x, T|b) over the hyper surface of hadronization 
determined by T(x, r) = Tc, and considering the decay of the excited states to the ground state, we can calculate 
the number A^aa of survived J/ips in a heavy ion collision. Let's first examine the differential nuclear modification 
factor RAA{pt) — NaaIpi)/ (-^c-^pp(pt)) as a function of transverse momentum pt in mid rapidity region, where 
Npp is the number of J/ips in a nucleon-nuclcon collision and Nc the number of nucleon-nucleon collisions. To see 
the largest mean field effect, we calculate first RAA{Pt) in central collisions with b = where the formed fireball 
is hot and large and the surviving time is long. FiglT] shows our numerical results for Au+Au collisions at top 
RHIC energy ^/snn = 200 GeV. The thick and thin lines indicate our results with and without considering the 

J/^ mass shift, by taking /j/^ and respectively. The dotted, dashed and solid lines are the calculations 

with initial production only, regeneration only and the total. 

From our numerical results, the series /j"^ [n = 0, 1, 2, • • •) converges rapidly, there is 

3% for any of the three cases. This confirms our qualitative estimation that the mean field effect is mainly reflected 
in the change in the threshold, the attractive force is a second order effect and the higher order corrections with 
n > 2 can be safely neglected. In all the numerical calculations in this paper we take /j/^ = f j/\p- Since the initial 
production has ceased before the formation of the hot medium, the change in the threshold does not affect it, and 
the correction is only from the small mean field force. That is the reason why the initial contribution shown in 
FiglHis almost the same for the two cases with and without considering the mass shift. However, the regeneration 
happens in the hot medium, it is affected by both the reduced threshold and the mean field force, the overall 
correction should be much larger in comparison with the correction to the initial production. Considering the fact 
that the regenerated J/^s in the early stage of the QGP will be eaten up by the hot medium and only the J/ips 
regenerated in the later stage of the QGP at lower temperature can survive, the enhanced regeneration is mainly 
in the low pt region, as shown in FigHJ At pt — 0, the total Raa goes up from 0.46 to 0.58, the enhancement 
is 26%! Since some J/ips are kicked to low pt region by the mean field force, there is a little reduction of Raa 
in the mid pt(2 — 3 GeV) region. At high enough pt, there is almost no effect of the mass shift, and the J/V' 
distribution is dominated by the perturbative QCD in the initial production. In comparison with our previous 
calculations 0, [sH] , the Gaussian smearing treatment [TtI . Isij for the Cronin effect used here reduces the Raa at 
high Pt. 
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FIG. 1: The differential nuclear modification factor RAA{pt) as a function of transverse momentum pt in Au+Au collisions 
at impact parameter b = and top RHIC energy ^snn ~ 200 GeV. The initial production, the regeneration, and the total 
are shown by dotted, dashed and solid lines. The thick and thin lines are the calculations with and without considering the 
mean field effect. 

To compare with the experimental data, we show in Figl2] the RHIC data at centrality 10% [s^, [s^l and our 
calculation at 6 = 4.2 fm. With decreasing centrality, the fireball temperature drops down and the mean field 
effect is gradually reduced. However, from 6 = to 4.2 fm, the change in the mean field is small, and the model 
calculations are almost the same. We show in Figl2]only the total calculation. Since the current data are still with 
large uncertainty, both the calculations with and without J/ip mass shift seem in reasonable agreement with the 
data. Especially, the dip at about pt — 2.5 GeV coming from the competition [sj] between the initial production 
and regeneration is slightly amplified by the attractive mean field force. To confirm the mean field effect, we need 
more precise data at low pt- At high pt where the mean field effect disappears, the higher dissociation temperature 
for fast moving charmonia (3^ will increase the J/tp surviving probability. 

We calculated also the RAA{Pt) in central Pb+Pb collisions at LHC energy ^/Inn = 2.76 TeV, as shown in 
FigEl Due to the higher temperature, larger size and longer life time of the formed fireball, most of the initially 
produced J/4>s at LHC are eaten up by the medium, leading to a very small contribution to the Raa at any pt- On 
the other hand, the regenerated J/ips in QGP at very high temperature T > Td are immediately dissociated in the 
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FIG. 2: The differential nuclear modification factor RAA{Pt) as a function of transverse momentum pt in Au+Au collisions 
at impact parameter b = 4.2 fm and top RHIC energy y'SNN ~ 200 GeV. The data are from [1^, [sj], and the thick and 
thin lines are the full calculations with and without considering the mean field effect. 



medium, and only the regeneration at temperature T < Td contributes to the finally observed J/^s. Therefore, 
the regeneration at LHC is not sensitive to the fireball temperature, and the mean field effect on J/^ production 
is similar to that at RHIC, see FiglS] 




FIG. 3: The differential nuclear modification factor RAA{pt) as a function of transverse momentum pt in Pb+Pb collisions 
at impact parameter 6 = and LHC energy -y/sjvjv = 2.76 TeV. The initial production, the regeneration, and the total are 
shown by dotted, dashed and solid lines, and the thick and thin lines are the calculations with and without considering the 
mean field effect. 



While the mean field effect changes significantly the differential nuclear modification factor at low transverse 
momentum, it does not affect the global yield remarkably. From Figs[T][31 the most important mean field effect is 
at = and the effective region is around pt — 0.5 GeV which is much less than the averaged J/tp transverse 
momentum (pt) ~ 2 — 3 GeV at RHIC and LHC energies. Therefore, the pt-integrated nuclear modification factor 
Raa{Np) as a function of the number Np of participant nucleons becomes not sensitive to the mean field effect. 
From our numerical calculations, the mass shift induced enhancement of Raa{Np) is very small in peripheral and 
mid-central collisions and less than 10% even in most central collisions. Different from Raa which is a summation 
of the initial production and regeneration, the averaged transverse momentum square {pD is governed by the ratio 
of regeneration to initial production [s^, [sE] . We found that the enhanced J/ ij) yield at low pt leads to a less than 
10% suppression of {pf) in most central collisions. 

In summery, we developed a perturbative expansion to study the mean field effect on J/?/; production in heavy 
ion collisions. Taking the mean field term as a perturbation, we analytically solved the J /ip transport equation 
with elastic and inelastic terms to any order and found a rapid convergence of the perturbative expansion. The 
initial J /ij} production, which happens before the formation of QGP, is not sensitive to the mean field force, but 
the continuous regeneration, which happens in the evolution of QGP, is significantly affected by the mean field. As 
a result, the differential nuclear modification factor of J/V' is strongly enhanced at low pt in heavy ion collisions 
at RHIC and LHC. 
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